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It is shown that the nearest neighbor coupling between the carbon atoms is not enough to repro- 
duce the phonon spectrum as observed in graphene, the second neighbor force constant is essential. 
For completeness we have rederived the phonon spectrum of one dimensional chains and two di- 
mensional lattices like square lattice and triangular lattice with nearest and next nearest neighbor 
coupling among the atoms. This article is essentially an account of the technical details of a review 
talk on phonon spectrum of graphene presented at Institute of Physics. 



I. INTRODUCTION 



Graphene is a one atom thick single layer of hexagonally arranged carbon atoms which has become experimentally 
accessible very recently. In this review only the in-plane vibration spectrum of perfect graphene has been calculated. 
The motivation was whether it is possible to find out the relevant features of phonon dispersion considering only nearest 
neighboring coupling among atoms because the important aspect of electronic dispersion in graphene is obtained with 
nearest neighbor hopping only, particularly the linear dispersion near the Brillouin zone corner (K point). Also in 
the electron-phonon coupling in this system only those electrons around the K point are going to take part. So it is 
relevant to see the nature of the vibration spectrum around the K point. In section(II) the preliminaries of lattice 
vibration within harmonic approximation have been discussed, then in section(III) we discuss about vibration spectra 
of standard ID monatomic linear chain, diatomic chain, 2D square lattice without basis and with basis, triangular 
lattice and finally the graphene. Section(IV) compares the results obtained here for graphene vibration with the 
results calculated by L. M. Woods and G. D. Mahan^ and L. A. Falkovsky^. 



II. SOME PRELIMINARIES 



The Hamiltonian of a vibrating crystal is given by 
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The notations used here are standard, that is, M n is the mass of the n th atom, Ui{n, I) is a small displacement of 
the n th atom in the I th cell along i th direction, i, j represent components of Cartesian coordinate axes and U is the 
ion-ion interaction potential. Moreover, <j>ij(m,n;l,l ) is defined as the force acting on the n th atom in the I th cell 
along i th direction due to a unit displacement of the m th atom in the V th cell along j th direction. The equation of 
motion of the n th atom in the I cell is given by 
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From the translational symmetry of a crystal the small displacement can be written as 
u i (n,l) = M- i u in e-<" t -*-^), 



(2) 



(3) 



where Ui„ is the amplitude of vibration along i direction of the n atom, u> is angular frequency, k is wave vector, 

~* —1/2 

R n (l) is the lattice translation vector and the factor M n has been chosen for convenience in further calculation. It 
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follows that 
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Dijimn, k)u mj , where D lJ (mn, k) = (M m M n )"5 ^ ^ ^'"^ (4) 



is the dynamical matrix which contains all the informations regarding the vibration of the lattice. Since this matrix is 
hermitian [Dij = -D^J, the eigenvalues are real and they give the vibration spectrum. Non-trivial solutions for these 
set of equations will be obtained if 

\Dij(mn, k) - Lo 2 SijS mn \ = 0. (5) 
We shall now discuss the properties of atomic force constant </>y(m, n;l,l ). 

1. Because </>jj(m, n; I, I ) is a second differential, 

2. Translational symmetry means that <f>ij(m, n;l,l ) depends only on R m (J) - Rn(l')- Hence, 

tpij ^ ^ — 4>ij q ' j ^ 7 where refers to an origin of coordinates. (7) 

3. The potential energy and the force on a given atom should be invariant under rigid body displacement of the 
whole crystal. This requires that 

E, fm,n\ „ . . 

n,Z x 7 

4. If the crystal has inversion symmetry then 

f m,n\ ( m, n\ , . 

These properties are useful in finding out relations between various force constants. To calculate the atomic force 
constants, let Uj(m) be the relative displacement of the m th atom in the j th direction with respect to the n th atom. 
Then the force acting on the n th atom in i th direction due to displacement of m th atom only is 

fi = 7mej(m) ej(m).Uj(m), where (10) 

3 

e(m) is the unit vector along R m (l) and 7 m the spring constant between n th and m th atoms. Then the total force 
acting on the n th atom in i th direction due to all other atoms is 



Fi = y^7 m ef(m)y^gj(m).Mj(m). (11) 

ml j 

comparing the above equation with equation^ and keeping in mind the definition of tpij (m, n; 0, 1), we see that 

--f m ei(m)ej{m). (12) 
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III. ONE AND TWO DIMENSIONAL LATTICES 
A. Monatomic chain 



We shall now consider the vibration of a monatomic linear chain shown in figure(l). Let 'a' be the lattice constant 
and 7 the nearest neighbor spring constant. Here the relevant unit vectors are 
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FIG. 1: Monatomic Chain. 



e\ = I, e 2 = -i. 

For one dimensional case the index (ij) corresponds to (xx) only and with one atom per cell the atomic index (m,n) 
in unimportant. So the indices in the force constant 0jj(m, n; 0, Z) may be reduced to 0(0, Z), where T represents 
nearest neighbor cell. According to the definition in equation (O, the force constants are calculated as 

0(0, 1) = -7 exx exx = -7 and 0(0,-1) = -7 e 2x e 2x = -7. 

Hence from the sum rule in equation (J5]), the self force constant is given by 

0(0,0) = -[0(0,1) + 0(0,-1)] = 2 7 . 

In this case the dynamical matrix is a (lxl) matrix and the clement is determined following the definition in equation 
O, i.e., 
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[0(0, 0) + 0(0, l)e lka + 0(0, - l)e- 4fea ] 
[1 — cos(fca)] . 

The eigen frequencies from equation ([3]) are given by 



w = 2J-^-sin(fca/2) 



(13) 



The brillouin zone of a one dimensional chain is a line with first brillouin zone boundaries at ±(n/a). The dispersion 
curve is shown in fig. (2). In the long wavelength limit the vibration spectrum starts from zero and goes linear; it 




FIG. 2: Vibration spectrum of a one dimensional linear chain, 
corresponds to the acoustic branch. 

B. Diatomic chain 

Let us consider a diatomic chain shown in fig. (3) with lattice constant 2a, where 'a' is interatomic distance. Mi 



FIG. 3: Diatomic Chain. 

and M 2 are the masses in the unit cell and 7 is the coupling constant between them. In this case the unit vectors are 
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same as that of monatomic chain. For linear chain the (ij) index may be omitted from 0y . From equation (|12p . the 
force constants are 



o, o) v \o,-i 



= -7- 



The symbol </>(l,2;0, 0) indicates force constant between first and second atoms in the th cell and ^(1,2; 0,-1) 
indicates the same between first and second atoms in the th and — 1 th cell. The self force constant <f>(l, 1; 0, 0) for 
the first atom is calculated from the sum rule in equation JHJ) as 
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Concentrating on the second atom in the unit cell, we get 
% 1\ _ /2, V 

v o, \) v \o, o y 

Hence from the sum rule 
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The dynamical matrix elements are 

D(22) . * and D(21) _ --^ (1 + „-*•) 
The eigcn frequencies from equation ([5]) are given by 



^=l(- + -)±J(- + -) - 4sin(fcaa) . (14) 
The dispersion curves are plotted graphically in fig. (4). Due to difference in masses between two atoms in a unit cell, 




FIG. 4: Phonon spectrum of a diatomic chain. 

there is a gap between the bands at the zone boundary ±(7r/2a) . If the two masses are equal, then the situation is 
basically one dimensional monatomic chain. So in that case there will be no gap at the boundary and the corresponding 
result will be in folded zone scheme, which is equivalent to the dispersion curve shown in figure (2) for monatomic 
chain. 
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C. MONATOMIC SQUARE LATTICE 

Now let us consider a two dimensional monatomic square lattice. The nearest neighbors are numbered as 1 to 4 
and next near neighbors as 5 to 8. The nearest neighbor and next nearest neighbor coupling constants arc 71 and 72 
respectively. The next step is to calculate a set of unit vectors for the different atoms and then the force constants 
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FIG. 5: Monatomic square lattice. 



can be determined. The unit vectors are 



-j, e 3 = -i, 
1 - 1 - A 1 , 1 - 

As we have numbered the neighboring atoms, it is convenient to write the force constants as <pij(m,n) instead of 
<t>j,j(m, n;0, I). According to equation (|12p . the force constants are 
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The self force constant for the atom is 

<j>xx(0, 0) = -[^(0, l) + (j)xx(0, 2) + <j> xx (0, 3) + 0^(0, 4) + <j) xx (0 : 5) + 
<f> xx (0, 6) + ra (O, 7) + xx (O, 8)] 
= 2 (71 + 72) . 

Similarly, 

^(0, 0) = 2(71+72). 
According to equation ([¥]) the dynamical matrix elements are 

D xx (k) = ^[^(0, l)e-^ a + ra (O, 3)e^ a + ^(0, 5 )e- , < fc « a+fc » a) + ^(0, 7 ) e -(-^-^) + 

^(0, 8)6-^-^+^+^(0, 6)e-^ a -^ a )+^(0, 0)] 
2 

= — [71(1 — cos k x a) +72(1 — cos k x a cos k y a)] . 



e\ = 1, e 2 ■ 
e*4 = j, e 5 
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Similarly, 

D xy (k) = D* x (k) = ^t[ — cos(k x a + k y a) + cos(k x a — k y a)] 

2 ^ 2 • U ■ U A 

= sin k. r a sin k„,a and 

M 
2 

Dyy(k) = — [71(1 — cos k y a) +72(1 — cos k x a cos k y a)}. 

The square brillouin zone of a square lattice with high symmetry points r(0, 0), X(n/a, 0) and L(ir/a, ir/a) is shown 
in fig. (6). The eigen frequencies from equation ([5]) considering atomic mass as unity along the symmetry directions 




FIG. 6: Brillouin Zone of square lattice. 



A = (T - X) arc 

loI = 2( 7 i +7 2 )(1 -cos k x a), (15) 

cl>2 = 272(1 — cos k x a); (16) 

along S = (T - L) 

Luf = 2 [71(1 — cosfc^a) + 272 sin 2 k x a\ , (17) 

uj 2 2 = 271(1 — cosfc x a); (18) 
and along Z = (X — L) 

w\ = 2 [2 7l + 72(1 + cosft B o)] J (19) 

uol = 2 [71(1 — coskyo) + 72(1 + cosfc^a)] . (20) 



The vibration spectra along various symmetry directions are shown in fig. (7). From the curve it is evident that phonon 




r x l r 



FIG. 7: Vibration spectrum of square lattice. 

frequencies are degenerate at T and L points whereas they are non-degenerate at X point. From the above equations 
it is evident that one of the phonon branches has zero frequencies along the T — X direction in absence of second 
neighbor coupling. But under same condition it has non-zero frequencies along X — L and T — L directions. 

D. SQUARE LATTICE WITH A BASIS 

Let us now consider a square lattice with a basis of two atoms of masses M\ and M 2 . The interatomic distance 
is 'a', i.e., lattice periodicity is 2a. In fig. (8) the two atoms in the th cell have been numbered as and 1. The 
atoms numbered as 1, 2, 3, 4 and 6, 7, 8, 9 are the nearest and next nearest neighbors respectively of the th atom. The 
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FIG. 8: Diatomic square lattice. 



nearest and next nearest neighbors of the first atom are 0, 5, 6, 7 and 3, 4, 11, 10 respectively. Let 7, 7! and 72 be the 
coupling constants among nearest Mi, M%; Mi, Mi and M2, M2 respectively. Proceeding in the same way like the 
monatomic square lattice, the different force constants become 

<M0, 1) = M (O, 2) = 4> yy {Q, 3) = 4> yy {Q, 4) = -7, 

<j) xx {0, 3) = M (O, 4) = cf> vv (0, 1) = TO (O, 4) = 0, 

<M0, 1) = 0^(0, 2) = <j> xy (0, 3) = 0^(0, 4) = 0, 

(t>xx(l, 0) = <f> xx (l, 5) = ytf (l, 6) = TO (1, 7) = -7, 

6) = 0^(1, 7) = W (1, 0) = (f)yy(l, 5) = 0. 

The corresponding (f> x s y and are zero. 

<j) xx {0, 6) = 0^(0, 7) = cf> xx (0, 8) = (f> xx {0, 9) = ~, 
ra (O, 6) = ra (O, 7) - ra (O, 8) = cf> yy (0, 9) = -y, 
<M0, 6) = 0^(0, 8) = 0^(0, 6) = 0^(0, 8) = - y, 
<^(0, 7) = <fe y (0, 9) = 0^(0, 7) = </v(0, 9) = y- 

Hence from the sum rule in equation ([5]) the self force constants for the th atom in the th unit cell are 

(f> xx (0, 0) = cf> yy (0, 0) = 2 7 + 2 7l and 
<M0, 0) = yx (O, 0) = 0. 

Concentrating on the second atom in the th cell, we get the following force constants 
<(> xx (l, 10) = <f> xx (l, 4) = <f> xx (l, 3) = XX (1, ll) = -y\ 

4>yy{l, 10) = 4>yy(l, 4) =0^(1, 3)=^J/(l, H) = ~, 

0^(1, 10) = XJ/ (1, 3)=^ x (l, 10) = <^ x (l, 3) = -y\ 
XJ/ (1, 11) = 0^(1, 4)=^ x (l, ll) = ^ x (l, 4) = y. 
The corresponding self force constants for the second atom are 
4>xx{l, 1) = TO (1, 1)^27 + 272 and 

<j> X y{l, 1) = <j>y x {l, 1) = 0. 

The dynamical matrix elements are 

£>xx(00) = ^[0 XX (O, 6) e -^ a+fc ^ + <M0, 8)e^ a+fe ^ + 

XX (O, 9)e l ^ a - fc « Q )+0 ra (O, 7)6-^-^+0^(0, 0)] 
2 

= — [7 + 71(1 — cosfc x acosfcj,a)] . 



8 



Similarly, 

271 271 

D xy (00) = -— sin k x a sin k y a, D yx (00) = -— sin k x a sin k y a, 

2 2 
D yy (00) = —[7 + 71(1 -cosfcajacos/fya)], D xx (ll) = — [7 + 72(1 - cos k x a cos ky a)], 
1V1\ 1VI2 

272 272 

D xy (ll) = -— sin k x a sin k y a, D yx (ll) = -— sin k x a sin k y a, 

2 27 

D yy( u ) = irr[7 + 72(l -cos k x a cos k y a)], D xx (01) = cos k x a, 

1VI2 v M1M2 

D xy (01) = D yx (01) = 0, D yy (01) = --^L=cosk y a, 

D xx (10) = — 7 =^=cosfc :c a ! D xy (10) = D yx (10) = and D yy (10) = — ^2= cos fc y a. 

y lVl\lVl2 \J M\Nl2 

The solutions to the secular equation ([5]) are plotted along the symmetry directions of the square brillouin zone with 
high symmetry points T(0, 0), X(0,ir/2a) and L(ir/2a, n/2a) for three different cases. 

Case 1 : Two masses in the basis are same and 7 7^ 0.0 but 71=72 = 0. In this case both acoustic and optical modes 




FIG. 9: Vibration spectrum of diatomic square lattice with equal masses in the basis. 



are degenerate at T point; at X point two are degenerate and two are nondegenerate but at L point all four are 
degenerate. Four branches are non-degenerate along T — X direction. In this direction one of the acoustic branches 
is having vanishing frequencies and one optical branch is dispcrsionless. Along X — L direction one acoustic and one 
optical branches are degenerate and along V — L direction both the branches arc degenerate. 

Case 2 : Two masses in the basis are same, 7 7^ 0.0 and 71 = 72 7^ 0.0. Here also both acoustic and optical modes 



M1=M2 




I X L 



FIG. 10: Vibration spectrum of diatomic square lattice with equal masses in the basis and with n. n. and n. n. n coupling. 

are degenerate at T point; at X and L points two are degenerate and two are nondegenerate. Four branches are 
non-degenerate along T — X direction as well as along X — L and L — T directions. 

Case 3 : Two masses in the basis are unequal and also 71 7^ 72 7^ 0.0. Here both acoustic and optical modes are 
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FIG. 11: Vibration spectrum of diatomic square lattice with unequal masses in the basis and with n. n. and n. n. n. coupling. 

degenerate at T point but at X and L points all are nondegenerate. All the four branches are nondegenerate along 
r — X, X — L and L — T directions. 

E. MONATOMIC TRIANGULAR LATTICE 

The structure of a triangular lattice is shown in the fig. (12). Each unit cell of this lattice consists of one atom. The 




FIG. 12: Triangular lattice. 



nearest neighbors of the atom are numbered as 1 to 6 and 7 is the spring constant among the neighbors. 
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are the unit vectors with respect to the atom. The force constants calculated by using equation p2[) are 
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Hence from the sum rule in equation |8]) the self force constants are <f> xx {Q, 0) = 37, 4> xy (Q, 0) = (f> yx (0, 0) = and 
4> yy (0, 0) = 37. The elements of the dynamical matrix are 



D xx {k) 
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k x a y^3k y a 

3—2 cos k x a — cos cos — 

2 2 



D xy (k) = D* yx {k) = - sin sin an d 
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Figure(13) shows the hexagonal brillouin zone of triangular lattice with high symmetry points T(0, 0), M(0, 2^/^30) 
and K (27r/3a, 2n/-\/3a). Considering atomic mass to be unity, the dispersion relations along the symmetry directions 
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FIG. 13: Brillouin Zone of Triangular lattice. 
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The graphical plot of the vibration spectrum is in fig. (14). Both the vibration spectra start from zero at the T point 




FIG. 14: Vibration Spectrum of Triangular lattice. 



(acoustic branch). They are nondegenerate at M point but degenerate at K point. The nearest neighbor coupling 
among atoms is sufficient to produce the vibration spectrum over the whole brillouin zone. 
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F. LATTICE VIBRATION OF HEXAGONAL GRAPHENE SHEET 

The structure of hexagonal graphene sheet is shown in fig. (15). It consists of two sublattices (say A and B) where 




FIG. 15: Structure of graphene sheet. 

A and B type of carbon atoms differ by their bond orientations. Thus each primitive cell contains two atoms. In 
the figure the two atoms in the th cell have been numbered as and 1. The nearest neighbors (n.n.) of one kind of 
atom (e.g. A) are three atoms belonging to other sublattice (B) and the next nearest neighbors (n.n.n.) are six atoms 
in the same sublattice. The atoms numbered as 1,2,3 and 4,5,6,7,8,9 are the nearest and next nearest neighbors 
respectively of the th atom. The nearest and next nearest neighbors of the first atom are 0, 4, 5 and 2, 3, 10, 11, 12, 13 
respectively. Let us now connect the n.n. and n.n.n. atoms by springs of spring constant 7 and 71 respectively. The 
distance between adjacent carbon atoms in the plane is ao(1.42A). The magnitude of the primitive vectors (ai and 
a 2 ) is 'a'(= aoV3). The set of unit vectors towards n.n. and n.n.n. with respect to the atom A (at 0) are the following 

1- \/3» A 1- V3- 
e l = 1, e 2 = --i + — j, e 3 = —-1 —j, 

n/3, 1- „ \/3~ 1 - 

n/3% 1~ 1- 
e7 = — 2~ l+ 2 h 68 = 2 J 69 = ~ J ' 

Similarly, the unit vectors with respect to the first atom can be calculated. According to the definition in equation 
(p~2|) . the force constants among nearest A — B and A — A type atoms are calculated. The force constants among 
nearest B — A and B — B type atoms are same as nearest A — B and A — A atoms respectively. 
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<M0, 7) = c(> xx (0, 8) = <j> xx (0, 4) = cf> xx {0, 5) = -^i, 

<M0, 7) = <^ x (0, 7) = cf> xy (0, 5) = ^(0, 5) = ^71, 

<^(0, 8) = ^(0, 8) = cj> xy {0, 4) = ^(0, 4) = -^7i, 
<M«, 7) = «^(0, 8) = cj> yy (0, 4) - ra (O, 5) - 

From sum rule in equation ([8|) the self force constants are 

<M0, 0) = -[^(0, l) + ^(0, 2) + 0^(0, 3) + ra (O, 6) + ^ x (0, 9) + 
<M«, 7)+<j> xx (Q, 8) + <t> xx (0, 5) + </> xx (0, 4)] 



^7 + 37i- 



Similarly, 



0^(0, 0) - <j) yx {0, 0) = and <j> yy {Q, 0) = -7 + 3 7l . 

For the second atom the force constants can be calculated in the same way and the self force constants are 

3 

<f>xx(l, l) = 0yt,(l) l) = ^7 + 37i and <f> xy (l, 1) = (f> yx (l, 1) = 0. 
Here (0 0) and (1 1) stand for A and B sublattices respectively. Now, the elements of the dynamical matrix are 
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3 . ,k v a 
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= D xx (ll) 
Let us call this element as A\ . Similarly, 

£^(00) = D yx (00) = D xy (ll) = D yx (ll) = - ll ^ s in(k x a^-) sin(^) = B l7 



D yy (00) = D yy (U) = -L 



M 



^7 + 7i{3 - cos(fc x a-^) cos(^) - 2 cos(k y a)} 



Dxx(01) = - 



_J_ 

M 

D* XX (W) 



<M0 l)e < ^ ) +4>xx(0 2)e l( ^l ^+^(0 3)e l( ^+^» 

Ci, 



e ~ vs ; + -e { ^> cos(-) 
2 v 2 ' 



V3 



D xy (01) = £,,(01) = D^(10) = D* yx (10) = -ij^ 7 e 1 ^ sin(^) = D x , 

D yy (01) = D^10)=-|^&co S (^)^ 2 . 

The secular determinant is 

A 1 - uj 2 B 1 d D x 

Bi A 2 -u? D x C 2 

C{ D\ Ai - uj 2 B x 

D\ CI B x A 2 - uj 2 

The hexagonal brillouin zone of graphene is shown in the fig. (16). In the first brillouin zone the symmetry 
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points arc r(0, 0), M(2n/a\/3, 0) and JsT(27r/av3, 2ir/3a). Let us now look for the eigen frequencies of lat- 
tice vibration at the above mentioned symmetry points of the brillouin zone within nearest neighbor coupling 
only, i.e. putting 71 = 0. Considering atomic mass to be unity, the eigenfrequencies at the T point are 
u)\ = 0, 1D2 = UJ3 = 0, UJ4 = \/3~7, at the M point are oj\ = y/j, uj-2 = \/2j, uj 3 = 0, 074 = y/Wy and at 

the K point are uj\ =0, w 2 — s/Sj/2, 0J3 = y / 37/2, = \/3^y. At Y point two of the eigenfrequencies are zero 
(acoustic branches). Other two non-zero values which correspond to optic branchesare are degenerate (which should 
be the case for a non-polar material). At M point all four branches are nondegenerate and at K point two of the 
branches are degenerate. The dispersion relations along symmetry line Y — M are 





3 7 

2 


- 7 


5 

4 " 


L (V ^ 

h cos(fc K a— 




3 7 
2 


+ 7 


5 

4 " 


L (i. ^' 
h cos(fc x a— 




= 


and 






2 
uJ 4 


- 3 7 . 









(26) 

(27) 

(28) 
(29) 



The dispersion curves along all the symmetry directions are shown graphically in fig. (17). Though in case of diatomic 




FIG. 17: Phonon dispersion in graphene with n.n coupling only. 

square lattice one of the acoustic branches vanishes only along r — X direction, in case of graphene one of the acoustic 
branches has totally zero frequencies over the entire symmetry directions when the next nearest neighbor coupling 
among the atoms is neglected. This indicates that to have quite satisfactory results over the entire brillouin zone, 
the n.n.n. coupling should be taken into account. In presence of both type of coupling, the eigen frequencies arc 
u>i =0, UJ2 = yf&Y, w 3 = 0, uj 4 — y/3j at the G point and wi = y/j + 671, w 2 = + 671, oj 3 = ^/tT, cj 4 = 
y/3j + 271 at the M point. The dispersion relations along all the symmetry lines in this case are shown in fig. (18). 



IV. RESULTS AND DISCUSSIONS 



Comparing the vibration spectrum of triangular lattice with that of graphene, we see that one optical branch and 
one acoustic branch are degenerate at the K point in graphene and both the acoustic branches are degenerate at K 
point in triangular lattice. We also see that the full phonon spectrum over the entire brillouin zone for triangular 
lattice is obtained with nearest neighbor coupling among the atoms but for graphene one of the acoustic branches is 
totally suppressed over the whole brillouin zone in absence of second neighbor interaction. In the work by Falkovsky^ 
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FIG. 18: Phonon dispersion in graphene along with n.n.n coupling. 

and earlier work by Woods and Mahan^ a detailed comparison with experimental result with graphene has been made. 
They describe the phonon dispersion with nearest and next nearest coupling among atoms. Woods and Mahani have 
considered a two parameter (viz. a and (3) model, where a is characteristic to two body central force between two 
neighboring atoms and it entirely depends on the separation between the atoms and (3 is characteristic to three body 
force related to bending of two adjacent bonds. Thus, (3 essentially gives the coupling between next-nearest neighbor 
atoms via the neighboring atom. Hence though carbon atoms in graphene arc bonded covalently, we have used an 
effective coupling between the second neighboring atoms (considering coupling constant very small compared to first 
neighbor coupling constant) instead of exactly following the microscopic details like them. It follows that at high 
symmetry points viz. T, M, K the analytical results are same in presence of nearest neighboring coupling only. In 
presence of two coupling constants we have got the analytic results at T and M points but not at K point. But 
the numerical plot shows that two optical branches are degenerate at T point, at M point all four branches are non- 
degenerate and at K point one acoustic and one optic branches are degenerate. These general features also match 
with the features obtained by Woods and Mahani. In Falkovsky's work^, unlike equation (fTS)) . the atomic force 
constants have been calculated by imposing constraints arising from the point group symmetry of graphene and the 
general features of vibration spectrum at the high symmetry points arc similar to results obtained here. Thus, we see 
that one acoustic branch of the vibration spectra of monatomic and diatomic square lattices is suppressed along one 
direction of the brillouin zone whereas in graphene it is suppressed over the entire brillouin zone in the absence of 
second neighbor coupling. So, to have a full dispersion curve over the whole brillouin zone in graphene as well as in 
square lattices, second neighbor coupling should be taken into account. 

ACKNOWLEDGMENT 

I am happy to thank Prof. S. G. Mishra for encouraging me and for his active help for this work, Dr. B. R. Sekhar 
and Mr. Nabyendu Das for fruitful discussions, Dr. Monodccp Chakraborty and Mr. Trilochan Bagarti for some 
technical helps. 



* Electronic address: rupali@iopb.res.in 

1 O. Madelung, Introduction to solid state theory. 

2 Landsberg, Solid State Theory Methods and Application. 

3 L.A. Falkovsky, |"arXiv:cond-mat 70702409^ 2. 

4 L. M. Woods and G. D. Mahan, Physical Review B, 61,10651(2000) 



